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ABSTRACT 

I demonstrate that precision timing of millisecond pulsars possess the capabilities of detecting the 
gravitational effects of intervening galactic substructure. This analysis is applicable to all types of 
collapsed baryons including stars, planets, and MACHOs, as well as many types of dark matter, in- 
cluding primordial black holes, scalar miniclusters, and sufficiently dense clumps of cold dark matter. 
The physical signal is quantified and decomposed into observable and unobservable components; tem- 
plates for the observable signals are also presented. Additionally, I calculate the expected changes in 
the observed period and period derivatives that will result from intervening matter. I find that pulsar 
timing is potentially a very useful tool for probing the nature of dark matter and to learn more about 
the substructure present within our galaxy. 
Subject headings: 



1. INTRODUCTION 

Our modern understanding of cosmology has revo- 
lutionized our picture of our own galaxy. What was 
once thought to be an island in the universe dominated 
by stars is now better modeled as a disk of mostly 
non-luminous baryons embedded in a much more mas- 
sive halo composed primarily of dark matter (' Cole et all 
[2005). However, the mass distribution and dumpiness of 
the both the disk and the halo, as well as the nature of 
the dark matter, are at present unknown. 

Baryons collapse to form n iolecular clouds, stars, an d 
planets, among other things (|Fukugita fc Peeblesl[200l . 
The location and abundance of the baryonic components 
of the galaxy are known primarily from their optical ef- 
fects, such as emission and absorption features. This 
tactic does not work for dark matter, though, since it 
does not inter act electroma^netic ally with any apprecia- 
ble strength (lAkerib et al.l l2005l ). However, both dark 
matter and baryons do interact gravitationally. Thus, if 
a galactic probe that is highly sensitive to the gravita- 
tional influence of the components of the galaxy can be 
identified, it will not only have the potential to detect 
the baryons that optical searches have missed thus far, 
but can also be used to detect dark matter within our 
own galaxy. 

Such a galactic probe highly sensitive to gravita- 
tional effects has been known sinc e 1982: the mil- 
lisecond pulsars (iBacker et al.l Il982l ). Over one hun- 
dred millisecond pulsars have already been identi- 
fied in our g;alaxy, bo th in the disk and the halo 
([Manchester et al.1 120051 ). Because of the stability and 
high accuracy of the timing measurements of millisecond 
pulsars, even small changes in the gravitational environ- 
ment as the light from a pulsar travels to Earth can, 
in princi ple, be detected (Larchenkova & Doroshenko 
fl995l:lFaTgiQn fc Conversa no 1997; Hos okawa et al.lll999l : 
i Siegel et al. |[2QQ7|; | Seto fc Cooray 2007 ) . 

With this paper, I seek to quantify the physical effects 
and observable signatures of baryonic and dark matter on 
the timing measurements of millisecond pulsars, with a 
view towards using millisecond pulsars to detect hereto- 



fore invisible baryons and dark matter. The layout of this 
paper is as follows: section 2 discusses the physical grav- 
itational effects of matter on the light-travel-time from 
a pulsar to Earth. Section 3 focuses on what compo- 
nent of the physical signal is uniquely ascribable to this 
gravitational effect, identifying points of confusion and 
their resolution. In section 4, templates for identifying 
transiting matter are created and presented, along with 
probabilities for detecting dark matter via this method. 
Finally, section 5 concludes this paper with a discussion 
of possible applications and future work as well as briefly 
addressing more complicated models of dark matter. 

2. GRAVITATIONAL EFFECTS OF MATTER 

When considering the gravitational effect of matter 
on the light-travel-time of a signal emitted from a pul- 
sar, one must consider the gravitational potential at ev- 
ery point along which the signal travels. A constant 
gravitational potential will result in no changes in the 
light-travel-time, but temporal variations in the grav- 
itational potential will lead t o time variatio ns due to 
the Shapiro time delay effect (^ Shapirol 11964 ). This ef- 
fect has been well-studied within our own solar system, 
as, for example, the moon causes periodic variations ev- 
ery twenty- eight days i n the arrival time of pulses from 
all pulsars ([Willi [19801 ) . It has also been used to dis- 
cover planets and binary companions to certain pulsars 
(Konacki & Wolszczan 2003). In this work, I address the 
possibility that there will be changes in the gravitational 
potential along the line-of-sight (LOS) due to intervening 
matter. Note that this phenomenon is physically simi- 
lar to the microlensing phenomenon, where a clump of 
matter causes a predictable change in brightness of a 
background source. Meaningful constraints have already 
been placed on the MAC HO content of the galactic halo 
from microlensing limits (|Alcock et al.lll998l ). However, I 
argue that pulsar timing has the potential to be an even 
more sensitive probe than microlensing of both collapsed 
baryons and, particularly, of non-baryonic dark matter. 

I consider the following configuration as illustrated in 
figure [U where a pulsar located at a fixed distance / from 
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Fig. 1. — Configuration for a transiting matter clump near the 
line-of-sight (LOS) from Earth to a pulsar located a distance I from 
Earth. Parallel components of position (x\\) and velocity (v\\) of 
the matter clump are along the LOS, while perpendicular compo- 
nents are in the plane perpendicular to the LOS. The perpendicular 
component of velocity is defined to be in the direction Vm, while 
the direction perpendicular to both the LOS and Vm is defined as 
the direction of the impact parameter, b. 

Earth has a clump close to the LOS of (presumably non- 
luminous) matter of mass Mdm- The position of the mas- 
sive clump is given by the components xy, which is the 
distance from Earth to the clump, and which is the 
perpendicular distance of from the LOS to the clump. xj_ 
can be further broken up into components {xm, where 
b is the impact parameter, and Xm is the component of 
position that will change over time due to the peculiar 
velocity v of the clump, which itself has components v\\ 
and Vm- 

Baryons make up approximately one-sixth of the mat- 
ter in the galaxy, wit h dark matter c omposing the other 
five-sixths (Wilkinso n fc Evanslll999l ). I estimate the to- 
tal mass of the galaxy tobe~1.2xl0^^ M©, and define 
/dm as the fraction of galactic matter (baryonic and dark) 
found in clumps of mass Mdm- Each clump is estimated 
to have a typical velocity of |iT| :^ 200kms~^ relative to 
the LOS to the pulsar, evenly distributed among its com- 
ponents, Vm and V|| . This is a reasonable estimate for the 
velocity relative to the LOS to the pulsar, as millisecond 
pulsars typically have large peculiar velocities relative to 
the bulk fiow of nearby baryonic matter. 

I assume henceforth that the clump of intervening mat- 
ter is either point-like or sufficiently far away from the 
LOS that the impact parameter (b) is greater than the 
physical radius of the clump of matter. Although this 
is not necessarily true for all forms of matter (WIMPy 
dark matter microhalos are the obvious exception), the 
calculations involved in computing the effects of extended 
dark matter microhalos are beyond the scope of this pa- 
per. With the clump of intervening matter treated as 
effectively point-like, the Shapiro time delay {At) due to 
any point-like mass Mdm is 
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where xy, x_\_^ and / are defined above. 

Equation ^ is the time delay induced by a clump of 
mass as compared with not having that mass present at 
all. However, physically, the mass is always present, and 
merely changes its position over time. In this scenario, 
the dark matter will transit from an initial position at 



an initial time (t^), and thus cause an initial time de- 
lay (Ati), to a final position (tf) with a final time delay 
(At/), leading to a physical delay (tdeiay) given in equa- 
tion 

^delay = Atf — Ati 
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If this physical delay were an entirely observable ef- 
fect, one would note immediately its linear dependence 
on mass (Mdm) and its logarithmic dependence on po- 
sition, x±. This would result in objects such as the 
Large Magellanic Cloud and Andromeda having much 
larger physical effects than anything transiting close to 
the LOS. The resolution to this problem, discussed fur- 
ther in section 3, is that not all of the physical effect 
is observable, and it is only the observable component 
that is of interest. I also note that the final expression in 
equation ([2]) is still completely analytically correct, but 
will prove useful when the special case of Vmtohs x± 
is examined below. 

3. PULSAR TIMING AND DETECTABILITY 

The change in the pulse arrival time, however, as given 
by equation ([2|), is the magnitude of the physical signal 
induced by a transiting clump of matter, not the magni- 
tude of the observable signal. The reason for this is that, 
a priori, the period and period derivative of the pulsar 
are not known. As such, when one measures the period 
(P) and the period derivative (P), one is not measuring 
the actual physical period and physical spin-up or spin- 
down of the pulsar. Instead, there are a number of effects 
that are confounded with the intrinsic period and period 
derivative, and they contribute, too, to what is observed 
for P and P. 

Any signal that contributes linearly with time to the 
arrival time of pulses will be confounded with the period, 
and any signal contributing quadratically with time will 
be confounded with the period derivative. The parallel 
component of the velocity of the pulsar relative to Earth 
is an example of an effect inseparable from the period. 
Because there is no way to measure the velocity of the 
pulsar independently, as the pulsar moves along the LOS, 
the distance along the LOS changes by an amount equal 
to the pulsar's velocity (vp) times the pulsar's period (P), 
causing a change in the apparent period of VpP/c^ where 
c is the speed of light. Likewise, the pulsar's parallel 
acceleration causes an analogous shift in the apparent 
period derivative, P. 

The way that pulsar astronomy accounts for these in- 
separable effects is to subtract them out. This is ac- 
complished through the measurement of the pulse arrival 
times of many pulses over considerably long temporal 
baselines, and constructing the best linear and quadratic 
fits to the pulse arrival times. From the linear fit and 
the uncertainty in it, the period (P) and the period un- 
certainty (AP) of the pulsar are determined. From the 
quadratic fit and its uncertainty, the period derivative 
(P) and its uncertainty (AP) are determined. 

What this means for transiting clumps of matter is 
that their contribution to the change in pulse arrival time 
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must be decomposed into a best fit linear + quadratic 
piece, which wih be absorbed into the period, period 
derivative, and their uncertainties. The remaining tim- 
ing residuals are the detectable portions of the physi- 
cal effect detailed in equation ([2]). It becomes immedi- 
ately clear that high-mass objects significantly far away 
from the LOS will have a significant impact on the ob- 
served period, but will have a negligible impact on the 
timing residuals. For objects significantly far away, i.e., 
xi_ ^ '^m^obs, it becomes useful to use the Mercator ex- 
pansion of the natural logarithm in equation ([2|) , obtain- 
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where the first term in the series will alter the period, the 
second will alter the period derivative, and only the third 
and higher terms will affect the timing residual. The 
effect on the timing residual, i.e., the detectable part, 
scales as {vm'tohs/xi_)^ ^ for objects far enough away that 
xi_ ^ '^m^obs is a valid approximation. As ('^m^obs/^±)^ 
is a number of order ~ 10~^^ for the Large Magellanic 
Cloud, it is therefore completely negligible. 

Since I have demonstrated that the effects of matter 
at too great a distance will be negligible in its effects on 
pulsar timing data, I am compelled to examine matter 
relatively close to the LOS. This begs the question "how 
close must one be?" to the LOS to have a shot at detect- 
ing a clump of mass Mdm- I approximate the answer by 
assuming the detectable component of the time delay to 
be the entire time delay of equation ([2]) with the first two 
terms of equation ([3]) subtracted out. I then numerically 
invert the equation to solve for the necessary (which 
will be equivalent to the necessary impact parameter, h) 
to be able to detect the matter clump in the timing resid- 
ual. It is assumed that the transiting clump of matter 
will pass through the point of closest approach during 
the timeframe of observation. (The point of closest ap- 
proach is also known as "conjunction" in some circles.) 
Figure [2] shows the maximum possible in parsecs that 
will result in a detectable signature in the timing residual 
as a function of mass (Mdm) of the matter and timing 
accuracy (tres) of the pulsar. 

Based on this derived distance, henceforth referred to 
as 6max, or the maximal impact parameter that will al- 
low for detection of transiting matter, a number of in- 
teresting quantities can be derived. Again, assuming the 
clump of matter of interest has a component of its ve- 
locity where Vm is in the plane perpendicular to the 
LOS, but is also orthogonal to the direction of the impact 
parameter, h. I assume that the pulsar is observed over 
a timescale tobs- In this case, there are four regions of 
interest to consider, as to whether the matter clump will 
pass within detectable range of the LOS, as illustrated in 
figure [3] below. For a clump originating in the innermost 
region containing point 1, the probability of detection is 
unity, whereas in the outermost region (containing point 
4), the probability is zero. In the region containing point 
2, at a radius 6max < r < Vmtohsi and the region con- 
taining point 3, at a radius Vm^obs < r < Vmtohs + ^max, 
the probabilities {p) are given by using the appropriate 
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Fig. 2. — The distance within which a clump of matter of mass 
Mdm must come of the hne-of-sight in order to be detectable in 
the timing residual. The five curves (left to right) are for pulsars 
that can deliver timing accuracies of 1 ns, 10 ns, 100 ns, 1 /iS, and 
10 /iS, respectively. This is for a relatively conservative estimate of 

Vmtohs — 10~^ PC. 

trigonometric laws in equation (|4|), assuming that the 
orientation of Vm is completely random. 
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Also note that equation (|4]) assumes that '^m^obs > ^max; 
if this is not true, then there is no region containing 
point 2, and the quoted probability for the region con- 
taining point three is valid for the entire range 6max < 

T <C Vrntohs H~ ^max* 

For each pulsar that is observed, the probability of 
observing a clump of matter passing within the critical 
distance to the LOS can be calculated. Every pulsar will 
have its own properties, including its own timing accu- 
racy (tres), period (P), period derivative (P), uncertain- 
ties in the period (AP) and period derivative (AP), and 
its distance from Earth (/). This probability, in addition 
to its dependence on tobs and will also depend on how 
long the pulsar is observed for (tobs), what the mass of 
the matter clump in question is (Mdm), the velocity of 
the transiting matter in the plane perpendicular to the 
LOS (vm)^ what fraction of the galaxy (by mass) is locked 
up in clumps of mass of Mdm (/dm), as well as the total 
mass of the galaxy, Mgai. 

The probability is calculated as follows: I as sume an 
NFW profile for the matter in our galaxy (Navar ro et al.l 
Il997l ). with a turnover radius at 25 kpc and our local po- 
sition located at 8 kpc, both measured from the galactic 
center. I further assume that a fraction /dm of the mat- 



Fig. 3. — Face on view of the LOS to a pulsar. Any clump of 
matter of mass M^m that comes within a distance 6max, as given by 
figure [21 will be detectable. The probability of a clump of matter 
coming within a distance 6max of the LOS for each of the four 
possible regions illustrated here are calculated in equation (|4]). 

ter is locked away in clumps of mass Mdm, which move 
randomly through the galaxy in all possible directions 
with a typical rms velocity of 200kms~^. This yields an 
rms velocity in the plane perpendicular to the LOS of 
141 kms~^. All volume in a cylindrical shape of length / 
and radius 6max is sensitive to transiting matter clumps. 

Figure [4] illustrates the probability of observing one or 
more clumps of transiting matter as a function of Mdm, 
tobs, and tres- is varied continuously, and the prob- 

ability {p) is shown for five different timing residuals of 
1 ns, 10 ns, 100 ns, 1 /is, and 10 jas are illustrated, along 
with four different observing times of 1, 10, 100, and 
1000 pulsar-years. / is fixed at 1 kpc, but the probability 
is directly proportional to I. The probability p is also 
directly proportional to /dm (assumed to be 1) and Mgai 
(assumed to be 1.2 x lO^^M©), and depends on Vm (as- 
sumed to be 141kms~^) in the same way as it depends 
on tobs (as shown in figured]). These probabilities are 
small and highly dependent on the above assumptions, 
but there may be ways to bring the total probability of 
detecting such matter close to unity, which I elaborate 
on in section 5. 

4. SIGNAL TEMPLATES AND FITTING 

More important than a calculation of the physical time 
delay induced by transiting matter, or even than the 
magnitude of the detectable component, it is vital to 
uncover what the signature of transiting matter will look 
like in the pulsar timing data. As the arriving pulses 
are very regular, it becomes very easy to calculate the 
physical signals in the timing residuals induced by the 
transiting matter, illustrated for a variety of impact pa- 
rameters in figure [5] for a clump of mass Mdm — 10~^Mq 
with a transverse velocity Vm = 141 kms~^, an initial po- 
sition = 0.002 pc, and an observing time (tobs) of 25 
years. 

As stated earlier, however, a signal like that shown in 
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Fig. 4. — Probability of detecting transiting matter of mass 
by looking at one given pulsar with timing accuracy tres for 
an observing time tobs- Short-dashed lines are for an observing 
time of 1000 years, long-dashed for 100 years, solid for 10 years, 
and dotted for 1 year. Lower masses can be probed with higher 
probabilities by improving the timing accuracy; the leftmost set of 
curves corresponds to a timing accuracy of 1 ns, and moving right, 
the following sets correspond to accuracies of 10 ns, 100 ns, 1/iS, 
and 10 /iS. Note that at sufficiently high masses, improvements in 
the timing accuracy offer no increase in probability. 

figure [5] will never show up as illustrated in the timing 
residuals. What is feasible, however, is to subtract off 
the best linear and quadratic fit, and extract the signal 
arising from that. One key point is that the best quadratic 
and linear fits to the data will change over time as more 
data is collected. Therefore, I have two possibilities at my 
disposal for searching for transiting matter with pulsars: 
one is to match the theoretical signal templates with the 
observed timing residuals, and the other is to look for 
systematic changes with observing time of the period and 
period derivative. I detail both methods here. 

I perform a least-squares fit for a second degree poly- 
nomial in tobs for the first halves of each of the curves 
(between 12 and 14 years of observing time) as well as 
for the full 25 years, calculating the linear and quadratic 
fits as well as the templates for the remaining signal that 
will be observable in the timing residual data. The three 
curves from figure [5] are analyzed as follows: the physi- 
cal time delay data is fitted to a second-degree polyno- 
mial using least-squares. The fit for the first half of the 
curve (up until the matter makes its closest approach 
to the LOS) is calculated first. Subsequently, the fit to 
the entire 25 year simulation is calculated. The differ- 
ences between the least-squares fit and the physical de- 
lay results in the signal that will show up in the timing 
residuals; all of these are illustrated in figure [6] below. It 
is immediately clear that observing the complete tran- 
sit, rather than just a partial component of it, results in 
a far greater magnitude of observable signal. Addition- 
ally, smaller impact parameters result in a greater signal 
strength. 
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Fig. 5. — The physical time delay induced by a transiting 
clump of matter of mass = 0-01 with an impact pa- 

rameter of 10~^ pc (thick solid line), 10~^ pc (thin solid line), or 
10~^pc (thin dashed line), velocity Vm = 141 kms"-*^, initial po- 
sition \x\ = 0.002 pc, and observed for 25 years. Clumps that 
pass even closer to the LOS will have greater physical signals, but 
are extraordinarily infrequent. Clumps that move with a greater 
Vm will have their timescales (along the x-axis) compressed, while 
those that are slower moving will have their timescales stretched. 
The signal strength is directly proportional to Mdm- 

As can also be seen from comparing the graphs on the 
left-hand side of figure [6] with those on the right-hand 
side, the least-squares fits change dramatically from the 
partial transit to the complete transit, as well as pro- 
ducing a greater observable signal in the residuals. How- 
ever, the differences in the polynomial fits correspond to 
a difference in the inferred period (P) and period deriva- 
tive (P). The differences in both P and P are compa- 
rable to the inherent uncertainties in P and P for many 
millisecond pulsars. For the case of an impact param- 
eter of 10~^ pc, I calculate a perceived change in P of 
4.1 X 10~^^s and a change in P of 3.8 x 10~^^, where I 
remind the reader that P is dimensionless. For an impact 
parameter of 10~^pc, P changes by 2.0 x 10~^^s and P 
changes by 5.3 x 10~^^. Finally, for an impact parameter 
of 10~^pc, P changes by 2.7 x 10~^^s and P changes 
by 7.0 X 10~^^. Changes in P and/or P over very long 
timescales of this magnitude may therefore be indicative 
of transiting clumps of matter. 

5. DISCUSSION 

I have demonstrated the quantifiable, observable ef- 
fects of transiting matter on pulsar timing measurements. 
Perhaps most importantly, I have crafted templates de- 
tailing the observable signal in the timing residuals as 
well as predicting the corresponding changes in P and 
P which accompany such a transit. These are effects 
which can lead to the indirect detection of either col- 
lapsed baryonic matter or dark matter in between Earth 
and a millisecond pulsar. 
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Fig. 6. — Templates for the signal that will show up in the pulsar 
timing residual data if there is a transiting clump of matter close 
to the LOS. The upper two figures are for an impact parameter 
of 10~^ pc, the middle two figures are for an impact parameter of 
10~^ pc, and the lower two figures illustrate an impact parameter of 
10~^ pc. The figures on the left demonstrate motion of the clump. 
The thin solid lines represent the physical time delay, the dotted 
lines represent the best second-order polynomial least-squares fit, 
and the thick solid line, the difference of the two, represents the 
signal that should be present in the pulsar timing data. This is for 
a transiting mass Mdm = 0-01 ^0 with velocity Vm = 141 kms"-*^; 
the magnitude of the signal scales directly proportional to mass, 
and the timescale along the x-axis scales proportionate to velocity. 

Current timing accuracies for millisecond pulsars are 
typically of order ~ 1 /is, although there do exist a num- 
ber of millisecond pulsars with accuracies of ^ 100 ns 
("KasDi et al.""l994"; 'van Straten et al.l l200ll : iJacobv etld] 
2005; Hotan et al. 2006), such as B1937-h21, J1909-3744, 
J0437-4715, and J1713+0747. Additionally, uncertain- 
ties in P are typically O{10~^^ s) and uncertainties in 
P are typically O(10~^^), although some are more accu- 
rate, for example, J1909-3744 has an uncertainty in P of 
2.0 X 10"^^ s and in P of 2.0 x 10"^^ ( Manchester etld] 
[2005h . While timing accuracies will have to improve 
by two orders of magnitude to detect clumps of mass 
^dm — 1O~^M0 with somewhat fortuitous transit pa- 
rameters, measurements of the uncertainties in P are 
nearly at the required level to give evidence of a transit 
today. (I remind the reader that the Mdm used in ear- 
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lier sections are 0.01 Mq.) I recommend that everyone 
working with mihisecond pulsars compare their timing 
residuals to these templates and also look for the corre- 
sponding systematic, coherent changes over time in the 
period and period derivatives of their pulsars. 

While the probability of observing a transit sufficiently 
close to the LOS in observing any one pulsar is exceed- 
ingly small, there are four reasons for optimism. First, 
the analysis presented is for one individual pulsar, but 
over 100 millisecond pulsars are presently known, with 
the total detectable galactic p ulsar population estimated 
to be approximately 30,000 (jKiel et al.ll2007f ). many of 
which will be discovered and timed to unp reced ented ac- 
curacy with the Square Kilometer Array (|Cord es 2005). 
Second, this assumes a distance to the pulsar of inter- 
est of 1 kpc. While this is a typical distance to the 
most accurate millisecond pulsars at present, millisecond 
pulsars have been detected as far away as the globular 
clusters, and new estimates are that millisecond pulsar 
timing is feasible for objects as far away as ~ 1.2 Mpc 
(|van Leeuwen fc Stap pers 2007), which would open up 
an extremely long baseline to millisecond pulsars in 
other galaxies (e.g., Andromeda), should this be correct. 
Third, the gravitational effects should be easily separable 
from other effects (e.g., dust lanes, hot gas, etc.) because 
of the greyness of the gravitational effects. While elec- 
tromagnetic effects caus e shifts of differing; mag:nitudes a s 
a function of frequency (|Weisberg fc Stanimirovi(^l2007l ) , 
the gravitational time delay is frequency independent. 
Other gravitational sources, such as gravitational scat- 
tering of the pulsar with a star or planet, or a gravita- 
tional wave effect on timing, will have a different tem- 
plate signal and can always be discriminated against by 
continued timing measurements. Finally, improvements 
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